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Q" 1 Abstract 

^jy The chaotical dynamics is studied in different Priedmann- Robertson- 

Walker cosmological models with scalar (inflaton) field and hydrody- 
namical matter. The topological entropy is calculated for some partic- 
^ ' ular cases. Suggested scheme can be easily generalized for wide class 

■ of models. Different methods of calculation of topological entropy are 

compared. 

I. INTRODUCTION 

In recent years the problem of chaos in general relativity and classical cos- 
mology has attracted great attention A special efforts has been made 
to resolve M the questions concerning the chaotical nature of the oscillatory 
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approach to singularity || or in other terms the chaotic nature of the mix- 
master universe [|J. However, already simple isotropical closed Friedmann- 
Robert son- Walker models manifest some elements of chaotical behavior which 
should be taken into account for a correct construction of quantum cosmo- 
logical theories |j| . 

The study of classical dynamics of closed isotropical cosmological model 
has a long history. First, it was noticed that in such models with a minimally 
coupled massive scalar field there is an opportunity to escape a singularity 
at contraction 0. Then the periodical trajectories escaping a singularity 
were studied || . In Ref . it was shown that the set of infinitely bouncing 
aperiodical trajectories has a fractal nature. Later this result in other terms 
was reproduced in our papers |I| |ll|, O . 



Here we would like to describe briefly the approach presented in [1C]. The 
main idea consisted of the fact that in the closed isotropical model with a 
minimally coupled massive scalar field all the trajectories have the point of 
maximal expansion. The localization of the points of maximal expansion 
on the configuration plane (a,0), where a is a cosmological radius, while (j) 
is a scalar field, could be found analytically. Then the trajectories could 
be classified according to localization of their points of maximal expansion. 
The area of the points of maximal expansion is located inside the so-called 
Euclidean or "classically forbidden" region. Numerical investigation shows 
that this area has a quasiperiodical structure, zones corresponding to the 
falling to singularity are intermingled with zones in which are placed points 
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of maximal expansion of trajectories having the so-called "bounce" or point of 
minimal contraction. Then studying the substructure of these zones from the 
point of view of the possibility of having two bounces, one can see that this 
substructure repeats on the qualitative level the structure of the whole region 
of the possible points of maximal expansion. Continuing this procedure ad 
infinitum one can see that as the result one has the fractal set of infinitely 
bouncing trajectories. 

The same scheme gives us an opportunity to see that there is also a set 
of periodical trajectories. All these periodical trajectories contains bounces 
intermingled with a series of oscillations of the value of the scalar field <fi. It is 
important that there are no restrictions on the lengths of series of oscillations 
in this case. In the paper of Cornish and Shellard || the topological entropy 
was calculated for this case by the methods of symbolic dynamics and it was 
shown that it is positive. Using symbolic dynamics is not new for models 
of theoretical physics. For example, in the well-known paper || it is used 
for investigation of the chaotic repellor, generated by the scattering of a 
point particle from three hard discs in a plane, which one can consider as a 
model of the time evolution of a metastable configuration of particles. This 
repellor consists of the trajectories which remain confined to the scattering 
region as t — > oo. If three hard discs with radius a are fixed in the plane 
at the vertices of an equilateral triangle with side R and R > 3a, then all 
the trajectories of this repellor are in one-to-one correspondence with doubly 
infinite sequences ...X-2X-1XQX1X2... of the symbols {1, 2, 3} with the natural 
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constraint x n+ \ 7^ x n . The symbolic dynamics for it is very simple: in the 
sequences only doublets such as 11, 22 and 33 do not occur. The topological 
entropy of the system equals In 2 and therefore such a repellor is chaotic. 



sible to generalize it for more complicated cases, but first we describe how 
introduction of matter or another modification of the model changes the 
structure of periodical trajectories. The structure of the paper is the follow- 
ing: in the second section we describe the properties of cosmological models 
under investigation, while in the third section we present the algorithm for 
the calculation of the topological entropy. 



First let us write down the action for the simplest cosmological model with 
the scalar field: 



where mp is the Planck mass, A is the cosmological constant. The equations 
of motion for a closed isotropic universe are 



Below we reproduce the calculations of Ref. || and show how it is pos- 



II. COSMOLOGICAL MODELS 




(2.1) 




(2.2) 
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+ m 2 = 0. 



(2.3) 
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The first integral of motion of our system is 




(2.4) 
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In the case when the cosmological constant is equal to zero, the form of the 
boundary of the Euclidean region is given by an equation which can be easily 
obtained from Eq. (|2.4| ): 

m V0 2 = — m 2 P . (2.5) 

47T 

The form of this Euclidean region is shown in Fig. 1(a). In this model, 
investigated in many papers || |], [10], [5| , there are periodical trajectories 



with an arbitrary number of oscillations of the scalar field. 

It can be understood by using two famous asymptotics for this dynamical 
system. For large the slow-roll regime leads to a quasiexponential growth 
of the scale factor (so-called inflation), and in the end of this regime 

a-end = a exp(2n((j) /M p ) 2 ) 

where do, (fio are the initial values and a en d is the scale factor in the end of 
inflation. 

For small the massive scalar field looks like dust perfect fluid but with 
small oscillations with the frequency equal to m. Using the Friedmann solu- 
tion for a closed dustlike universe, it is possible to express the scale factor in 
the point of maximal expansion a me through the value of the scale factor at 
the end of inflation a en d- 

(ma end ) 3 ~ ma me 
and the duration of the Friedmann stage 
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So, starting with O , 0o >> m P, the universe will go through 

N = m 3 a 3 end = m 3 aoexp(67r(0 o /M p ) 2 ) 

oscillations. This value can be arbitrarily large, because the value of the 
scalar field at the Euclidean boundary can also be arbitrarily large. 

Of course, it in necessary to prove that the bounce is indeed possible for 
any large value a me . It was done analytically by Starobinsky in [0 with the 
estimation of the bounce probability depending on a me . 

Inclusion of the positive cosmological constant investigated in JTI] implies 



two opportunities: if A is small in comparison with the mass square of the 
scalar field m 2 , the qualitative behavior is the same as in the model with- 
out cosmological constant; if the cosmological constant is of order m 2 the 
chaotical dynamics disappears in a jumplike manner [nj]. 

It is interesting from a calculational point of view to consider the case of 
the negative cosmological constant. In this case the Euclidean region has the 
form presented in Fig. 1(b) and numerical calculation shows that the possible 
number of oscillations is restricted. With the growth of the absolute value of 
the cosmological constant this number is decreasing and at |A| ~ 0.34m 2 all 
the periodical trajectories disappear and dynamics become regular. 

Let us include into consideration hydrodynamical matter with the equa- 
tion of state 

V = 7e, (2.6) 

where p is the pressure, e is an energy density, and 7 is the constant (7 = 
corresponds to dust matter, 7=1/3 corresponds to radiation, while 7 = 1 
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describes the massless scalar field). Our analysis is valid for 7 > —1/3. 
In this case the form of the boundary of the Euclidean region is given by 
equation 

m V0 2 = -^-m 2 P - — , (2.7) 

where .D is the constant characterizing the quantity of the given type of mat- 
ter in the universe and q = 3(7 + 1) — 2. The form of the Euclidean region 
described by Eq. (|2.7j ) is shown in Fig. 1(c). One can see that in this case 
the Euclidean region is restricted from above. Numerical calculation shows 
that in this case again only a restricted number of oscillations is possible. 
Moreover, the structure of periodical trajectories in this case is much more 
complicated. Indeed, the law is the following long series of oscillations that 
can come after the bounce, followed by a short series of oscillations, while 
the behavior of the trajectory after a short series of oscillations is less re- 
stricted. The concrete laws governing the structure of trajectories depend on 
the parameters of the model under consideration. However, in this case also 
it is possible to calculate a topological entropy, which will be demonstrated 
below. Completing a discussion about the model with hydrodynamical mat- 
ter one has to say that a large amount of matter suppresses chaotic behavior 
and at some critical values of the constant D periodical trajectories disap- 
pear and all the trajectories go to the singularity. For the above-mentioned 
cases, interesting from a physical point of view, values are the following: for 
the case q = 4 (massless scalar field), one has 

Dm 4 > 0.0028m 2 ,; 
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for q = 2 (radiation), 



Dm 2 > 0.0093m 2 



Pi 



and for q = 1 (dust matter), 

Dm > 0.023m 2 P . 

It is necessary to underline that bounces are still possible in such a sit- 
uation (unless D is much bigger than a critical value), but there are no 
periodical trajectories, because immediately after bounce the trajectory has 
a point of maximal expansion and restores their travel to singularity 

Now we would like to consider the last example: the cosmological model 



with a complex scalar field which was studied earlier in papers [fT3" |. Using 
the most natural representation of the complex scalar field 

<f> = xexp(i9), (2.8) 

where x is an absolute value of complex scalar field while 9 is its phase. 
This phase is cyclical variable corresponding to the conserved quantity - the 
classical charge of the universe, which plays the role of the quasifundamental 
constant of the theory: 

Q = a 3 x 2 6. (2.9) 
The shape of the Euclidean region presented in Fig. 1(d) is given now by 

2 2 2 3 2 Q 2 /r> i n\ 

m a x = — m P — -. 2.10 

47r a A x 2 

Again, as in the preceding example at some critical value of the charge Q 
periodical trajectories and chaotical dynamics vanish, while at lower values of 



charge there are again complicated rules governing the structure of periodical 
trajectories. The condition of regular behavior is given by the inequality: 

Qm 2 > 0.056m^. 

It is important to notice that the real structure of trajectories in the case 
of isotropic models with scalar field and matter or with complex scalar field 
is more complicated because there are two types of bounces. These bounces 
could be called "upper" and "lower" depending on the value of the scalar field 
at which they take place. The value of the scalar field separating the two 
types of bounces correspond to a periodic trajectory with a maximal number 
of oscillations of the scalar field. However, in what follows we shall not 
discriminate between the two types. Generalization of suggested algorithms 
for the calculation of topological entropy for the "alphabet" including two 
types of bounces is straightforward. 

III. CALCULATION OF THE TOPOLOGI- 
CAL ENTROPY 

Now, let us calculate the topological entropy for some of the considered cases. 
In our situation strange repellor Q is the set of all perpetually bouncing uni- 
verses. We shall find topological entropy of the Poincare return map S with 
discrete time (see Ref. 01), generated by the surface of section X consist- 
ing of the points, in which bounce or oscillation occurs: if x is the point of 
intersection (crossing) of the trajectory with the surface X, then Sx is the 
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a 



Figure 1: In Fig. 1 the shapes of Euclidean regions for different isotropic 
cosmological models with the scalar field are presented; a) simplest model 
with the scalar field, b) model with negative cosmological constant, c) model 
with hydrodynamical matter, d) model with complex scalar field. 
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point of its following intersection with X (by definition of Q each trajec- 
tory of Q crosses X infinitely many times). Such an entropy is independent 
of phase-space coordinates. We note that the topological entropy depends 
on the choice of Poincare surface for the section. Our Poincare surface of 
section X is natural for cosmological models with a massive scalar field for 
closed Friedmann-Robertson- Walker universes and plays the same role as the 
boundary of discs or obstacles for billiard systems. 

Topological entropy Ht "measures" the chaos of the dynamical system. 
There are many equivalent definitions of topological entropy (via finite par- 



tition of the phase space) in mathematics (see Ref. []15|). But for many 
dynamical systems (in particular, for the finite topological Markov chain, to 
which the symbolic systems considered below belong) this quantity coincides 
with the exponent of growth in the number of periodic points as their period 
increases. We remark that the set of periodic points and the set of all closed 
orbits are everywhere dense in space X and in the repellor Q, respectively, 
and thus, they contain basic information about dynamics. 

More precisely, let (S, X) be a dynamical system with discrete time. If 
x G X and S k x = x, the point x is called a periodic point of period k. We 
denote by N(k) the number of periodic points of period or length k. Then, 
the definition of topological entropy is as follows: 

H T = Urn sup ^ In N(k). (3.1) 

If Ht > 0, one can conclude that the dynamics is chaotic. 

One can quantify the length of the orbit by the number of symbols. To 
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begin with let us reproduce in some detail calculations from Ref. ||. It was 
a discrete coding of orbits, using two symbols: 

A, which is the bounce of the trajectory, and 

B, which is the crossing of line = 0. 

It is important that for our repellor for a given admissible infinite sequence 
of symbols A and B there is only one trajectory, which goes through bounces 
and oscillations in the order defined by this sequence. Therefore, in this 
sequel we can deal with only such sequences. 

For the simplest model with the scalar field there is the only prohibition 
rule: two letters A cannot stay together, which means that it is impossible 
to have two bounces one after another without oscillations between them 
(this condition is analogous to the impossibility of having two successive 
collisions about the same disc for the scattering from three discs). Thus our 
symbolic system is a topological Markov chain (of the first order) with two 
letters (see corresponding definition in Ref. |16|). Each topological Markov 
chain of the first order with k letters {B\, ...,Bk} is defined by topological 
transition matrix D. It is a k x k matrix of zeros and ones, the set of indices 
of which coincides with the set of letters. The entry Db 1 b 2 °f t ne matrix D 
is zero precisely when BiB 2 is a prohibited word of length 2. In our case, 
the topological transition matrix D is 




(the zero in this matrix corresponds to the prohibition of the word AA). 
There are many practical means of counting the topological entropy. For 
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example, calculating the trace of the matrix D k , we obtain the number of all 
period k points of the map S. Since Tr D k = X k + A§, where Ai and A2 are 
eigenvalues of the matrix D, it is sufficient to find the spectrum of D. In fact, 
we must calculate only the largest positive eigenvalue, which exists by the 
Perron- Frobenius theorem (see Ref . [17] ; such an eigenvalue is called Perron- 



Frobenius's number). In our case Ai = A^ 1 = ( v / 5 + l)/2. The logarithm of 
this number is sought to define the topological entropy Another means of 
calculating the topological entropy consists of finding the smallest positive 
root of the equation <Pa(z) — 1, where ^Pa( z ) is a generating function of 
periodic points without intermediate recurrences in the initial symbol A (see 
Ref. [IB |) In our case the next periodic points are of the following type: ABA, 



ABBA, ABBBBB...BBA, Thus tp A (z) = z 2 +z 3 + . .. + ... = z 2 /(l-z). 

The smallest positive root of the equation z 2 /(l — z) = 1 is (VE- l)/2. The 
logarithm of the reciprocal from this number is again a topological entropy. 

Now we present the method of recurrent relations, which as one can see 
in sequel is more successful for our purposes. Let us denote by Q(k) the 
number of "words" (trajectories) of length k satisfying this rule, which begin 
with A and end with A, and by P(k) the number of words which begins with 
A and end with B. (We remark that our definition of coefficients Q(k) and 
P(k) is distinguished from that of Ref. |J: in their paper P(k) is the number 
of period k points ending in A, and Q(k) is the number of period k points 
ending in B). Then one can easily write down the recurrent relations: 

Q{k + l)=P{k), 
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P(k + l) = Q(k) + P{k). (3.2) 

It is easy to get from Eqs. ( |3.2|) the following relation for P(k): 

P(k + l) = P(k) + P(k-l). (3.3) 

One can easily calculate that 

P(2) = 1, P(3) = 1 (3.4) 

and that Eq. ( |3.3j ) defines the series of Fibonacci numbers. Let us recall how 
to find the formula for the general term of the Fibonacci series. One can 
look for P(k) as a linear combination of terms X k , where A is the solution of 
equation 

X k+i = x k + X k-i 
or, equivalently, because we are interested only in nonzero roots 

A 2 - A - 1 = 0. (3.5) 

Looking for P(k) in the form 

P(k) = Cl \ k 1+ c 2 \ k 2 , (3.6) 

where Ai and A 2 are the roots of Eq. (|3.5|) and satisfying the conditions (|3.4j ) 
one can get 

PW = it 



+ (-1 



,k-2 



V5 



(3.7) 
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Substituting Eq.( |3.7| ) in the definition of the topological entropy ( |3.1| ) one 
can find 

H T = In f > 0, (3.8) 

where ( is the famous golden mean. It is clear that only the largest 
root of Eq. (|3.5|) is essential for the calculation of the topological entropy. 

We would like to point out that the number N(k) of all period k points 
grows as exp^Hrk) for k — > oo. In our case the number Q(k + 1) is in 
fact the number of all real closed orbits (with respect to dynamical system 
with continuous time, generated by Einstein's equations), which are coded 
by k symbols (here we distinguishes closed orbit and it's multiple). In other 
words, Q(k + 1) is the number of all periodic orbits of period k (each periodic 
point x of period k generates periodic orbit {x, Sx, S k ~ 1 x}\ it is clear 
that points Sx, S k ~ 1 x generate the same periodic orbit; thus the number 
of all period k points is much greater than that of all periodic orbits of 
period k) and its multiples. This number grows as /j,({xq = A})exp(HTk) 
for k — > oo, where [i is the measure of maximal entropy for our topological 
Markov chain and {xq = A} is the set of all sequences, which start with the 
letter A (so-called "one-dimensional" cylindrical set) |I9| . Such a measure is 
characterized by the next important property: it is the most chaotic measure 
among all invariant probability measures for our topological Markov chain 
(more precisely, it is a measure for which the Kolmogorov-Sinai entropy is 
equal to the topological entropy). For example, for the set of all infinite 
sequences of the letters A and B, such a measure is the Bernoulli measure v 
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with v({xo = A}) = u({x = B}) = 1/2. In the case of a topological Markov 
chain it is a Markov measure. A measure with maximal entropy plays an 
important role in the distribution of periodic points. There are many means 
of counting the probability of "one (and multi)-dimensional" cylindrical sets 
for such a measure. In particular, one can find the positive eigenvectors £ and 
£* for Perron- Frobenius's eigenvalue of the topological transition matrix A 
and its conjugate A*, respectively. Taking the product of the corresponding 
components of the vectors £ and £* and then normalizing it, we obtain the 
probability of each "one-dimensional" cylindrical set. On the other hand, 
this number is coefficient before the k power of the largest eigenvalue in 
the decomposition of Q(k + 1) = P(k). Thus, in our case it is equal to 

ci = (VE-i)/2VZ. 

Now one can go to a more involved case of the cosmological model with 
scalar field and negative cosmological constant. As was described above, in 
this model the periodical trajectories can have only a restricted number of 
oscillations of the scalar field 0. This rule can be encoded in the prohibition 
to have more that n letters B staying together, where the number n depends 
on the parameters of this model (this dependence obtained by numerical 
calculations is represented in Fig. 2). 

Thus we have a topological Markov chain of order n + 1. Correspondingly, 
the order of the topological transition matrix D{n) is equal to 2 n+1 . We 
see that the first method is not well applicable to our case (in general, to 
any typical big matrix; this method is useful for theoretical investigations). 
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Figure 2: In Fig. 2 the dependence of the number of possible oscillations TV 
on the ratio R = \A\/m 2 is represented. 
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n-1 

Using the second method we have <pa{z) = z 2+k — z 2 (z n — l)/(z — 1) 

fc=0 

and therefore topological entropy is equal to — In z(n), where z(n) is the 
smallest root of the next equation z n+2 — z 2 — z + 1 = (the largest positive 
eigenvalue A(n) of the matrix D(n) is connected with the number z(n) by 
the next relation: X(n) = ;?(n) _1 ). 

At this time the recurrent relations is as follows: 

Q{k + l)=P(k), 

P(k + 1) = Q(k) + P(k) - Q(k - n)9(k - n), (3.9) 

where the 9 - function is defined in the usual manner. We are interested in 
the limit k — > oo and can substitute instead of 0(k — n) number 1. Now one 
can write down the recurrent relation: 

P(k + 1) = P(k) + P(k -l)-P(k-n- 1), (3.10) 

which in turn implies the following equation for topological entropy: 

\{n) n+2 - X(n) n+1 - \{n) n + 1 = 0, (3.11) 

where topological entropy is equal to the logarithm of the biggest root of Eq. 

H T — In A(n). (3.12) 



For the small values of n, the biggest root of Eq. ( 3.11 ) could be found 
analytically. 

For n — 1, we have A(l) = 1 and topological entropy is equal to zero and 
chaotic behavior is absent, which is clear from physical point of view. 



For n = 2, we obtain 




1/3 (|(9 + V69)) 



1/2 



32/3 



1.32. 



(3.13) 



For n = 3, one can get 



a(3) = ! + !(—_ 

K ' 3 3 V 2 2 J 

1 /29 3 V / 93\ 1/2 , 
+- — + — — « 1.47. 

3 v 2 2 ; 



1/3 



(3.14) 



For higher values of n one can find A numerically, for example for n = 4, 



For large values of n one can find an asymptotical value for the biggest 
root A: 



Thus A(n) converges to golden mean exponentially fast when n — » 00 (or, 
negative cosmological constant tends to zero). 

Concerning the coefficient before /c-power of the largest eigenvalue or, 
the probability of "one- dimensional" cylindrical set corresponding to the 
letter A we may show that in our case it is equal to the reciprocal of 



z£(f A (z)\ z=z{n) = 2 + - n f; M (see Ref. [p|). Since z(n) tends 



the convergence of /i n ({x = A}) to yu({x = A}), where /x n is a measure of 
maximal entropy for the topological Markov chain with the topological tran- 
sition matrix D(n). Thus two basic parameters in the distribution of periodic 



A(4) ss 1.53. 




(3.15) 
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orbits of our approximating systems develop very regular behavior when the 
negative cosmological constant tends to zero. 

As it has already been described above in the model with the scalar field 
and matter or in the model with the complex scalar field and nonzero classical 
charge, the rules governing the structure are rather complicated. Neverthe- 
less, in this case also it is possible to calculate the topological entropy. Here 
we shall consider one particular example, however, the algorithm which shall 
be presented could be used for different sets of rules as well. 

Thus, let us formulate rules for our model: (1) It is impossible to have 
more than 19 letters B; (2) after a series with 19 letters B (and letter A), one 
can have the next series only with 1 letter B; (3) after a series with 18 letters 
B, one can have the next series with 1 or 2 letters B; (4) after a series with 
17 letters B, one can have the next series with 1, 2 or 3 letters B, . . . (20) 
after a series with 1 letter B, one can have the series with n letters B, where 
< n < 19. 

These rules define the topological Markov chain of order 22. The method 
of calculating the topological entropy via finding the spectrum of the topo- 
logical transition matrix and the method of generating functions are not 
applicable to this case. But the method of recurrent relations allows us to 
write the equation for topological entropy (we do not calculate the coefficient 
before the k power of the largest eigenvalue in the distribution of periodic 
orbits because it is a very difficult problem now). 

Let us notice that the system of rules has a remarkable symmetry with re- 



20 



spect to the number nc = 10, which makes further calculations more simple. 
These symmetrical rules give us a good approximation for the description 
of a real physical situation. The value ric is apparently a function of the 
parameters of the model under investigation. More detailed numerical inves- 
tigation implies a more complicated system of rules, however, it also could 
be formalized in the system of recurrent relations. Below we shall see that 
already the symmetric system of rules and a comparatively small number, 
which we have chosen for our example {nc = 10) gives rather a cumbersome 
equation for the topological entropy. 

Let us now introduce the following notation: 
Q(k) is the number of words which begin with letter A and end with letter 
A, 

Qi(k) is the number of words beginning with letter A and ending with a 
series with 1 letter B, 

Q^{k) is the number of words beginning with letter A and ending with a 
series with 2 letters B, 
. . . , 

Qig(k) is the number of words beginning with letter A and ending with a 
series with 19 letters B. 

The system of recurrence relations for these quantities is as follows: 

Q(k + 1) = Q 1 (k) + Q 2 (k) + ■■■ + Q 19 (k), 
Q 1 (k) = Q(k-l), 
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Qw(k) = Qi(fc-20), 

Qd(fc) = Q d -i(A; - 1) - Q2i- d (k - d - 1), 2 < d < 10 

Qd(fc) = g rf+ i(A; + 1) + Q 20 -d(fc - d - 1), 11 < d < 18 (3.16) 

Resolving this system with respect to Q(k), one get the following recurrent 
relation: 

Q(k + 1) = Q(k - 1) + Q(fc - 2) + Q(fc - 3) + Q(fc - 4) 
+Q(k - 5) + Q(fc - 6) + Q(k - 7) + Q(fc - 8) + Q(k - 9) 
+Q(A; - 10) + Q(k - 13) + 2Q(k - 14) + 3Q(k - 15) + 4Q(A; - 16) 
+5Q(k - 17) + 6Q(k - 18) + 7Q(k - 19) + 8g(A; - 20) + 9Q(k - 21) 
-9Q(k - 24) - 16Q(k - 25) - 21Q(ib - 26) - 24Q(A; - 27) 
-25Q(A; - 28) - 24Q(A; - 29) - 21Q(k - 30) - l6Q(k - 31) 
-9Q(k - 32) - 8Q(k - 36) - 21Q(k - 37) - 36Q(A; - 38) 
-50Q(A; - 39) - 60Q(A; - 40) - 63g(A; - 41) - 56Q(A; - 42) 
-36Q(A; - 43) + 36Q(A; - 47) + 8AQ(k - 48) + 126Q(A; - 49) 
+150Q(A; - 50) + 150QO - 51) + 126g(A; - 52) + 84Q(A; - 53) 
+36Q(A; - 54) + 28Q(A; - 59) + 84Q(A; - 60) + 150Q(A; - 61) 
+200Q(A; - 62) + 210Q(A; - 63) + 168g(A; - 64) + 84Q(A; - 65) 
-84Q(A; - 70) - 224Q(A; - 71) - 350Q(A; - 72) - 400Q(A; - 73) 
-350Q(A; - 74) - 224Q(A; - 75) - 84g(A; - 76) - 56Q(A; - 82) 
-175Q(Jfe - 83) - 300g(A; - 84) - 350g(A; - 85) - 280Q(A; - 86) 
-126QO - 87) + 126g(A; - 93) + 350g(A; - 94) + 525Q(A; - 95) 
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+525Q(A; - 96) + 350Q(A; - 97) + 126Q(fc - 98) + 70Q(k - 105) 
+210Q(A; - 106) + 315Q(A; - 107) + 280Q(A; - 108) + 126Q(A; - 109) 
-126Q(k - 116) - 336Q(A; - 117) - UlQ(k - 118) - 336Q(A; - 119) 
-126Q(k - 120) - 56Q(k - 128) - 147Q(A; - 129) - 168Q(A; - 130) 
-SAQ(k - 131) + UQ(k - 139) + 196Q(A; - 140) + 196Q(A; - 141) 
+84Q(A; - 142) + 28Q(A; - 151) + 56Q(A; - 152) + 36Q(A; - 153) 
-36Q(A; - 162) - 64Q(A; - 163) - 36Q(A; - 164) - 8Q(k - 174) 
-9Q(k - 175) + 9Q(k - 185) + 9Q(k - 186) + Q(k - 197) 
-Q(Jfc- 208), (3.17) 

which in turn gives the following equation for the topological entropy: 

209 _ 207 _ 206 _ 205 _ 204 _ 203 _ 202 _ 201 _ 200 

Ju tlj tlj «X «X Ju Ju tlj Ju 

■x 199 - x 198 - x 195 - 2x 194 - 3z 193 - 4a; 192 - 5a- 1 " 1 - 6x 190 

_ 7a; 189 _ gx 188 _ 9x 187 + 9x 184 + ^183 + ^182 + ^181 

+25x 180 + 24x 179 + 21a; 178 + lQx 177 + 9x 176 + 8x 172 + 21a; 171 
+36a; 170 + 50a; 169 + 60a; 168 + 63a; 167 + 56a; 166 + 36a; 165 - 36a; 161 
-84a; 160 - 126a; 159 - 150a; 158 - 150a; 157 - 126a; 156 - 84a; 155 
-36a; 154 - 28a; 149 - 84a; 148 - 150a; 147 - 200a; 146 - 210a; 145 
-168a; 144 - 84a; 143 + 84a; 142 + 224a; 137 + 350a; 136 + 400a; 135 
+350a; 134 + 224a; 133 + 84a; 132 + 56a; 126 + 175a; 125 + 300a; 124 
+350x 123 + 280a; 122 + 126a; 121 - 126a; 115 - 350a; 114 - 525a; 113 
-525a; 112 - 350a; 111 - 126a; 110 - 70a; 103 - 210a; 102 - 315a; 101 
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+ 336x' 



.89 



+ 1262 s8 + 56x 8U + 147x 79 + 168x 78 + 84x 77 - 84a;' 



-196x bS - 196x b ' - 84x bb - 28x 5 ' - 56x 5b - 36x 



,55 



+36a; 4b + 64a; 45 + 36a; 44 + 8x d4 + 9x dd - 9x 
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(3.18) 



Resolving numerically Eq. (|3.18| ) one can find the biggest root which is equal 



to 



A 



1.61771. 



Correspondingly the topological entropy is given by the logarithm of the 
biggest root. We remark that the degree (209) of our equation is much 
smaller than the order of the topological transition matrix (2 22 ). It means 
that most of the eigenvalues of the topological transition matrix are equal 
to zero and, therefore, this matrix contains much useless information about 
our dynamics. Thus we see that the method of recurrent relations is more 
economic in this case. 

One can easily see that similar calculations can be done for every value of 
ric- Here we can give the list of numbers A corresponding to different values 

of Tic' 



n c = 2; A « 1.37747 



n c = 3; A ~ 1. 



51714 



n c = 4; A w 1.57388 



n c = 5; 



A 



1.59837 
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n c = 10; A m 1.61771 

n c = oo; A = (1 + V5)/2 « 1.61803. 

Thus, the topological entropy tends to the golden mean rather rapidly. 

Using this algorithm it is possible to get the equations for the topological 
entropy in the case of the two types of bounce (see the end of Sec. 2). We 
will not present here the explicit forms of these equations, which differ from 
the previous case only in the values of nonzero coefficients, but present only 
some results. 

First, the analog of the golden mean will be the number 2. We can see 
that there is not a smooth transition in the limit D — > (when D = we 
have the simplest case with only one type of bounce). But from a physical 
point of view, when D — > the values of the scalar field, corresponding to the 
upper bounce grow, cross the Planck limit, and the upper bounce becomes 
impossible. So the transition to the simplest case takes place. 

And finally, we present for illustration some values for the topological 
entropy in the symmetric case with two types of bounce: 



n c = 2; 


A R 


i 1.82657 


n c = 3; 


A R 


i 1.94561 


n c = 4; 


A R 


i 1.98273 


n c = 10; 


A f 


a 1.99999 



nc = oo; A = 2. 
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We investigate the behavior of the topological entropy of some cosmolog- 
ical models, when they are close to the simple isotropical closed Friedmann- 
Robert son- Walker model [its parameters (negative cosmological constant, 
charge) tend to zero]. It was shown that chaotic behavior is robust and 
topological entropy converges to that of the FRW model. In the case of the 
presence of a negative cosmological constant, we give a full description of the 
corresponding symbolic dynamics and the laws of the convergence of basic 
parameters in the distribution of periodic orbits to that of the FRW model. 
The situation of the presence of matter in the model with a real scalar field 
and the model with a complex scalar field are more difficult, but using the 
method of recurrent relations we may find the values of the topological en- 
tropy. This is a case when the topological entropy is calculated for such com- 
plicated symbolic dynamics (the corresponding topological transition matrix 
has 209 nonzero eigenvalues). The scheme described in this example can 
be applied to many different physical models, obeying the different sets of 
prohibition rules governing the structure of periodical trajectories. 
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